Thermal Conductivity of Nanotubes Revisited: Effects of Chirality, Isotope Impurity, 

Tube Length, and Temperature 
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We study the dependence of thermal conductivity of single walled nanotubes (SWNT) on chirality, 
isotope impurity, tube length and temperature by nonequilibrium molecular dynamics method with 
accurate potentials. It is found that, contrary to electronic conductivity, the thermal conductivity 
is insensitive to the chirality. The isotope impurity, however, can reduce the thermal conductiv- 
ity up to 60% and change the temperature dependence behavior. We also found that the tube 
length dependence of thermal conductivity is different for nanotubes of different radius at different 
temperatures. 

PACS numbers: 65.80+n, 66.70.+I, 44.10.+i 



Carbon nanotube is one of exciting nano-scale mate- 
rials discovered in the last decade. It reveals many ex- 
cellent mechanical, thermal and electronic properties^. 
Depends on its chirality 2 ** 3 -, the nanotube can be either 
metallic or semiconducting. For example, for zigzag (9,0) 
and (10,0) tubes, their radius are almost the same, but 
(9,0) tube behaves metallic and (10,0) tube semiconduct- 
ing. At room temperature, the electronic resistivity is 
about 10~ 4 — 10 -3 ficm for the metallic nanotubes, while 
the resistivity is about 10f2cm for semiconducting tubes 4 . 
The 10% difference in radius induces the change of elec- 
tronic conductivity in four orders of magnitude. One may 
ask, whether thermal conductivity is also very sensitive 
to the chirality like its electronic counterpart? 

On the other hand, the isotope impurity reduces ther- 
mal conductivity of most materials, such as germanium 
and diamond*^. It is surprising that 1% 13 C in diamond 
leads to a reduction of thermal conductivity up to 30% 7 . 
Is there same effect in carbon nanotubes? 

Moreover, in electronic conductance, SWNT reveals 
many ID characters 4 . However in thermal conduction, 
it is still not clear whether the conduction behavior is 
like that one of ID lattice or a quasi ID (ID lattice with 
transverse motions) or that one in a 2D lattice. 

These questions and many other relevant proper- 
ties of nanotubes are very important and should be 
understood before the nanotubes are put into any 
practical application. Indeed, recent years have wit- 
nessed increasing interesting in thermal conductivity 
of nanotubeo 8 ! 9 ! 10 ! 11 ! 12 ! 13 . 14 ! 15 ! 16 ! 17 ! 18 ! 19 , the questions 
raised here are still open. 

In this paper, we study the effects of the chiral- 
ity, isotope impurity, tube length and temperature 
on SWNTs' thermal conductivity by using the non- 
equilibrium molecular dynamics (MD) method with bond 
order potential. This approach is valid as it shows that 
at finite temperature, phonon has a dominating contribu- 
tion to thermal conduction than electron does 20,21 . We 



should pointed out that the thermal conductivity calcu- 
lated in this paper is exclusively from lattice vibration. 
Of course, for the metallic nanotubes, the electrons may 
give some but limited contributions 20 * 2 ^, but this is not 
the main concern in our paper. 

The Hamiltonian of the carbon SWNT is: 



H=y 

^ \ 2m, 



V u 



(1) 



where V io = fc(nj)[V R (nj) + byVAfaj)] is tne Tersoff 
empirical bond order potential. Vn(rij), and V^r^ ) are 
the repulsive and attractive parts of the potential, and 
f c {r) depending on the distance between atoms, by are 
the so-called bond parameters depending on the bound- 
ing environment around atoms i and j, they implicitly 
contain many-body information. Tersoff potential has 
been used to study thermal properties of carbon nan- 
otubes successfully 1 ^* 2 ^. For detailed information please 
refer to Ref2£. 

In order to establish a temperature gradient, the two 
end layers of nanotube are put into contact with two 
Nose-Hoover heat bathes* 2 ^ with temperature Tl and Tr 
for the left end and the right end, respectively. Free 
boundary condition is used. All results given in this pa- 
per are obtained by averaging about 10 6 ~ 10 7 femtosec- 
ond (fs) after a sufficient long transient time (usually 
10 6 ~ 10 7 fs) when a non-equilibrium stationary state is 
set up and the heat flux, J, becomes a constant. The 
thermal conductivity, k, is calculated from the Fourier 
law, 



J 



-kVT, 



(2) 



where J is defined as the energy transported along the 
tube in unit time through unit cross section area, and 
VT = dT/dx is the temperature gradient. In this paper, 
we chose d — 1.444 as the tube thickness, thus the cross 
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FIG. 1: The temperature profiles of (9,0) and (10,0) SWNT at 
300A" . The temperature is obtained by averaging over 5 x 10' 
fs after dropping the 10 6 - 1 n7 
length, L, is 1081. 
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10 fs transient time. The tube 



TABLE I: Thermal conductivity for different tubes. 



Tube 


(9,0) 


(10,0) 


(5,5) 


K (W/mK) 


880 


770 


960 


Radius, r (A) 


3.57 


3.97 


3.43 


Relative thermal conductance 


1.0 


0.97 


1.05 



section is 2nrd, where r is radius of the tube. In the fol- 
lowing we shall discuss the chirality dependence, isotope 
impurity and tube length effect. 

Chirality Dependence. The thermal conductivity of 
zigzag and armchair SWNTs of same length but with 
different radius are calculated. Fig. ^ shows the tem- 
perature profiles of (9,0) and (10,0) nanotubes at 300A". 
These two temperature profiles are very close to each 
other, thus the temperature gradient, dT/dx, is almost 
the same. The thermal conductivity, tube radius and 
relative thermal conductance are listed in Table 1. The 
difference in thermal conductivity comes mainly from ra- 
dius difference. If we use thermal conductance (defined 
as thermal conductivity times cross section area), the 
relative value (to (9,0) tube) for (9,0), (10,0) and (5,5) 
SWNTs with the same length is 1 : 0.97 : 1.05. 

It is clear that, unlike its electronic counterpart, the 
thermal conductivity/conductance of SWNTs does not 
depend on the chirality and/or atomic geometry sen- 
sitively both at low temperature and room tempera- 
ture. Our MD results are consistent with that one from 
Landauer transmission theory 21 . The electron DOS of 
SWNTs depends on chirality. There is an energy gap at 
Fermi level in (10,0) tube, while no such a gap in (9,0) 
tube. However, the phonon DOS's in different tubes do 
not show any significant difference 2 ^. 

Isotope impurity effect. Isotope impurity affects many 
physical properties of materials, such as thermal, elastic, 
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FIG. 2: (a) Thermal conductivity, ft, versus 14 C impurity 
percentage for a (5,5) SWNT at 300A". (b) Thermal conduc- 
tivity, k, versus temperature for a (5,5) pure 12 C nanotube ( 
solid A) and a (5,5) SWNT with 40% 14 C impurity (•). The 
curves are drawn to guide the eye. 



and vibrational properties 2 ^. There are three isotopes 
of carbon element, 12 C, 13 C, and 14 C. They have the 
same electronic structure, but different masses. Here we 
study the effect of 14 C impurity on thermal conductivity 
of SWNTs. In our calculations, 14 C atoms are randomly 
distributed in a 12 C SWNT. To suppress the possible 
fluctuations arising from random distribution, an average 
over 10 realizations is performed for each conductivity 
calculation. We also do the calculation for tube with 13 C 
impurity and similar effect is found. 

Fig. 13a) shows the dependence of thermal conductiv- 
ity on the impurity percentage. We select armchair (5, 5) 
tube with a fixed length of about 60A The thermal 
conductivity decreases as the percentage of 14 C impurity 
increases. With 40% - 50% 14 C, the thermal conductiv- 
ity is reduced to about 40% of that one in a pure 12 C 
SWNT. This is similar to the isotope effect on thermal 
conductivity of diamond 7 . 

The thermal conductivity decreases more quickly at 
low percentage range than at high range. From this 
curve, we can estimate roughly that the thermal con- 
ductivity decreases about 20% with only 5% 14 C isotope 
impurity. This decrease is not as rapid as that one in a 
diamond that 1% isotope impurity can reduce thermal 
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conductivity as much as 30%t 

This result tells us that one can modulate the thermal 
conductivity of carbon nanotubes by adding 14 C or other 
isotope impurity as it alters only the thermal conductiv- 
ity and has no effect on the electronic properties. 

The thermal conductivity, k, for (5,5) pure 12 C nan- 
otube and that one with 40% 14 C impurity at different 
temperatures with same tube length are shown in Fig. 
Hfb). The difference for these two cases are obvious. The 
isotope impurity changes completely the temperature de- 
pendence behavior of thermal conductivity. For a pure 
tube (solid A), there is a maximum at about Tm ~ 250K. 
Below this temperature, n increases when T is increased. 
Above Tm, k decreases with increasing T. However, in 
the case with isotope impurity, there is no maximum in 
the curve. The thermal conductivity monotonically de- 
creases as the temperature increases. 

These phenomena can be understood from the phonon 
scattering mechanism. Increasing temperature has two 
effects on thermal conductivity. On the one hand, the 
increase of temperature will excite more high frequency 
phonons that enhance thermal conductivity. We call this 
effect "positive" effect. On the other hand, the increase 
of temperature will also increase phonon-phonon scatter- 
ing that in turn will increase the thermal resistance, thus 
suppress the energy transfer. We call it 'negative" ef- 
fect. The thermal conductivity is determined by these 
two effects that compete with each other. 

For a pure SWNT, at low temperature regime, the 
phonon density is small, the "positive" effect dominates, 
thus the thermal conductivity increases with increasing 
temperature. However, at high temperature regime, as 
more and more (high frequency) phonons are exited, the 
"negative" effect dominates, which results in the decrease 
of thermal conductivity as temperature is increased. The 
results in Fig[3fb) is consistent the results from Savas et 
alA 

However in the case with impurity, the scattering 
mechanism changes. In this case, most high frequency 
phonons are localized due to the impurity. The main 
contribution to heat conduction comes from the low en- 
ergy phonon that has long wavelength. The "positive" 
effect is largely suppressed, and in the whole tempera- 
ture regime, the "negative" effect dominates that leads 
to a decrease of thermal conductivity as the temperature 
is increased, as is seen in Fig^b). 

Tube length effect. Recent years' study on heat conduc- 
tion in low dimensional lattices shown that for a one di- 
mensional lattice without on-site potential, thus momen- 
tum is conserved, the thermal conductivity k diverges 
with system size (length)2£, L, as k ~ L 13 , with (3 = 2/5. 
If the transverse motions are allowed, like in the quasi 
ID case2£, then (3 — 1/3. Moreover, it has been found 
that this anomalous conduction is connected with the 
anomalous diffusion^. For more detailed discussion on 
anomalous heat conduction and anomalous diffusion in 
different models such as lattice models, billiard gas chan- 
nels, and nanotubes, see Ref^S. However, in a 2D system, 
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FIG. 3: The temperature profile of (5,5) SWNT and ID lat- 
tice at (a) 2K, (b) 300K, and (c) 800K. The temperature 
is obtained by averaging over 5 x 10 6 fs after dropping the 
10 6 ~ 10 7 fs transient time. 

it is still not known at all (both analytically and numeri- 
cally) that what shall be the divergence form: power law 
form or logarithmic form or anything else. 

As for the SWNT, Maruyamai^ has studied the k{L) 
for tubes of different radius, and found thati^ the value 
of (3 decreases from 0.27 for (5,5) tube, to 0.15 for (8,8) 
and 0.11 for (10,10). 

Here we investigate this problem from different as- 
pects, namely, we study the change of (3 in different tem- 
peratures. For comparison, we also study a ID carbon 
lattice model with the same interatomic potential. The 
carbon-carbon atom distance is also 1.44A 

The temperature profiles for a carbon SWNT and a ID 
carbon lattice are shown in Fig [3] for different tempera- 
tures. In this figure, both SWNT and ID lattice has 100 
layers of atoms. FigOJa) demonstrates clearly that at low 
temperature as low as 2K, there is no temperature gradi- 
ent in both the SWNT and ID lattice. This resembles the 
ID lattice model with a harmonic interaction potential 31 . 
This can be understood from the Taylor expansion of the 
Tersoff potential by keeping up to the second order term. 
Because at low temperature, the vibrations of atoms are 
very small, the potential can be approximated by a har- 
monic one. And in SWNT, the vibration displacement in 
transverse direction is much smaller than the one along 
the tube axis and can be negligible. This result means 



4 



40 

30 
20 

10 



fa) 


A, 




(5,5) 4 




(10,10) 


300 K > 




300 K / 


/A 






: 4/ L 040 




/ L 0,6 


,,l 





30 



20 



10 



^ 7 
I 6 

5 





(d) , 


(5,5) f 


(10,10) ^ 


800 K / 


800 K 






/ L ° 26 


/ T 0.12 . 


L 





1 10 100 

Tube length (nm) 

FIG. 4: The thermal conductivity, k, versus tube length, L, 
in log-log scale for (5,5) and (10,10) tubes at 300K and 800K. 
In all cases, k ~ L 13 with (3 changes from case to case. The 
solid line, whose slope is the value of (3, is the best fit one. 



that energy transports ballistically in SWNTs at low tem- 
perature. 

However at high temperature, the situation changes. 
As is shown in Fig. I^b) for 300K, and Fig. EJc) for 
800K, there is still no temperature gradient in the ID 
lattice, but in the SWNT, temperature gradient is set up. 
In ID lattice with Tersoff potential, the increase of tem- 
perature does not change the harmonic character; while 
in the SWNT at room temperature and higher temper- 
ature, transverse vibration increases as temperature in- 
creases. The temperature gradient is established due to 
the interaction of the transverse modes and longitudinal 
modes. 

The thermal conductivity, k, versus the tube length, L, 
is shown in log-log scale in Fig^Ja-d) for (5,5) and (10,10) 
SWNT at 300K and 800K, respectively. Obviously, the 
value of P depends on temperature as well as tube radius. 
For a SWNT, (3 decreases as temperature increases (cf. 

(a) and (c), (b) and (d)); and at the same temperature, 
(3 decreases as the tube radius increases (cf. (a) and 

(b) , (c) and (d)). This can be qualitatively explained 



by the modes coupling theory 28 . At high temperature, 
the transverse vibrations are much larger than that at 
low temperature, thus interaction between the transverse 
modes and longitudinal modes becomes stronger, which 
leads to a smaller value of (3. For (5,5) tube at T=300 K, 
because the small tube radius, the transverse modes can 
help set up the temperature gradient but the coupling 
with the longitudinal mode is still very weak, and the 
longitudinal modes dominate the heat conduction, this 
is why in this case the thermal conduction behaviour is 
very close to that one in the ID Fcrmi-Pasta-Ulam type 
lattices 2 ^, namely, the thermal conductivity, k, diverges 
with L, as L 0A . 

It is worth pointing out that the absolute values of 
thermal conductivities given in Fig. 4 for (5,5) and 
(10,10) tubes are different from that ones given in 
Refsii^i^. The reasons are that: (a) the wall thickness 
used in Refsii^i^ was 3.4 A, while in our case we use 
1.44 A, this leads to our results are at least about 2.5 
times larger than those ones in Refs 1 ^* 1 ^; (b) the ther- 
malstate we used are different from that one in Refsi 1 ^ 1 ^; 
(c) boundary conditions are different, we use free bound- 
ary condition while Refs 1 ^^ use fixed boundary condi- 
tion. If we use conductance (get rid of the effect of tube 
thickness and tube radius) rather than the conductivity, 
we believe that our results are consistent with that one 
from others 1 ^ 1 ^* 1 ^. Our value of divergent exponent for 
(5,5) SWNT at 300K is slight different from Maruyama's 
on o 14 i 15 . This might be caused by different boundary 
conditions and different heat baths used. 



In summary, we have studied the effects of chirality, 
isotope impurity, tube length, and tube radius on thermal 
conductivity in SWNTs. Our results show that the ther- 
mal conductivity is insensitive to the chirality. However, 
the introduction of isotope impurity suppresses thermal 
conductivity of carbon nanotubes up to 60% and change 
the temperature dependence behavior. Moreover, at low 
temperature the heat energy transfers ballistically like 
that one in a ID harmonic lattice, while at high temper- 
ature, thermal conductivity diverges with tube length, 
L, as L 13 . The value of (3 depends on temperature and 
tube radius. This unique structural characteristic makes 
SWNT an ideal candidate for testing heat conduction 
theory, in particular, the mode-coupling theory 2 ^ in low- 
dimensional systems. 
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